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Abstract 



The Concurrent Track Evolution method, which was introduced in a previous 
paper 0, has been further explored by applying it to the propagation of track 
candidates into an inhomogeneous magnetic field volume equipped with track- 
ing detectors, as is typical for forward B spectrometers like HERA-B or LHCb. 
Compared to the field-free case, the method was extended to three-dimensional 
propagation, with special measures necessary to achieve fast transport in the pres- 
ence of a fast-varying magnetic field. The performance of the method is tested 
on HERA-B Monte Carlo events with full detector simulation and a realistic 
spectrometer geometry 
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1 Introduction 



Dedicated B physics experiments for hadronic interactions, as HERA-B [fj] and 
LHCb |3j are designed as forward spectrometers because of the kinematics of the 
b hadrons, whose decay particles are produced with a large Lorentz boost. Even 
with highly granular tracking devices, as honeycomb drift chambers and strip 
detectors, the large track density leads to high cell occupancies, which can range 
up to 20% in the "hottest" parts of the HERA-B outer tracker. Under these 
circumstances, pattern recognition becomes a crucial issue. 

In a previous paper Jjj], the Concurrent Track Evolution strategy has been 
presented as a method of track finding in a multi-planar tracker geometry. Con- 
current Track Evolution is a track following method based on the Kalman fil- 
ter ||, [|, H [7J], which evaluates the available paths for a track candidate con- 
currently to find the optimal solution, but also keeps the combinatorics at a 
reasonable level (see H for details). The method has been successfully applied to 
find initial track segments in the geometry of the HERA-B pattern tracker, which 
consists of four superlayers of inner and outer tracking devices in the field-free 
region between the magnet and the RICH. The layout of the HERA-B spectrom- 
eter is shown in fig. [l|. A more detailed explanation of the HERA-B tracker and 
the reconstruction strategy can be found in [|TJ. 

The next step in the HERA-B track reconstruction is the track propagation 
through the magnetic field. Typical sizes of the magnetic field components B x , 
By and B z as a function of z are displayed in fig. 0. The main bending component 
(By) has a bell-shaped dependence on z, and there is clearly no significant subset 
of the magnet tracking system in which the field can be regarded as homogeneous. 
Moreover, the x and z components are sizeable, making it impossible to restrict 
pattern recognition to projections while making full use of the track model. The 
necessity to deal with five track parameters simultaneously, and the additional 
hits generated by spiralling particles make standalone track finding within the 
magnet a problematic task. Instead, the HERA-B track reconstruction concept 
operates by first finding straight-line track segments in the pattern tracker (de- 
scribed in [|J) and then propagating them upstream through the magnet area. In 
a subsequent step, these track candidates will have to be matched to the track 
segments in the silicon vertex detector. A considerable fraction of the Kg from 
the golden decay mode B° — > J/ipK s decays too late to produce a sufficient num- 
ber of hits in the acceptance of the vertex detector, so that the vertex must be 
reconstructed from the main tracker information. 

Compared to pattern recognition from scratchin the magnet, this propagation 
approach has the advantage that only one track parameter, the momentum, is 
poorly known and must be fitted during the propagation. There are however 
several challenges involved: 

• Due to the non-linear nature of the fit, iteration, in general, may be required 
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to obtain the optimal result. In particular, it might happen that upstream 
hits are not found within the predicted error margin because a wrong initial 
momentum value in the downstream section caused derivatives to be taken 
at the wrong positions in parameter space. Iteration of the various track 
candidates in the context of Concurrent Track Evolution, on the other hand, 
would multiply the computational effort to a worrisome degree. 

• The inhomogeneity of the field necessitates a much more costly, numerical 
transport of the track parameters than in the field-free case. The Kalman 
filter requires also the transport of the covariance matrices, for which the 
derivative matrices of the transported versus the original parameters must 
be calculated. 

• The different superlayer structures of outer and inner magnet tracker rep- 
resent additional challenges to the navigation concept. 

2 Geometry 

As an example, we will use the current HERA-B detector layout as it is presently 
implemented in the geometry definition of the experiment. This layout is shown in 
fig. [I]. On the right-hand side, four superlayers of inner and outer tracking devices 
form the pattern tracker, which is used for finding the initial straight-line track 
segments. The left part shows the magnet tracker through which these segments 
have to be propagated. The magnet tracker consists of eight superlayers of outer 
(MC01. . .MC08) and four superlayers of inner tracker (MS01, MS03, MS05, MS07). 
The outer tracker is based on honeycomb drift chamber technology with 5 and 
10 mm diameters, with an assumed resolution of 200/xm (see |l|, [| for further 
details). Each outer tracker superlayer in the magnet consists of three tracking 
layers with the orientations +a, 0°, —a with a stereo angle of a ~ 100 mrad, 
except for MC01 which has four layers in a 0°, —a, 0°, +a arrangement to satisfy 
requirements of second level triggering, and MC05 which consists only of a single 
0° layer because of space restrictions inside the magnet yoke. The inner tracker is 
built of micro-strip gaseous chambers (MSGC) with a typical resolution of 80/iin. 
Each superlayer consists of one zero degree and one stereo layer with a stereo 
angle of +100 or —100 mrad, alternatingly, except for MS01 which has four layers 
with the sequence 0°, —a, +a, 0°. The sector structure of the magnet tracker 
superlayers is similar to that of the pattern tracker (see 

Compared to the setup shown in the stations MC09 and MS09 have been 
removed because they are already in an area of weak magnetic field and are 
not expected to improve the track parameters already measured by the pattern 
tracker significantly^. The gaps in the second half of the magnet tracker are 

4 It is also foreseen to remove MC07 
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reserved for pad chambers designed to generate a pretrigger for high momentum 
tracks, which is intended to select B° — > 7r + 7r~ decays. In front of MS01, an 
additional silicon strip detector superlayer is foreseen to cover the innermost 
angular region. Both the pad chambers and the silicon superlayer were not used 
in this study. 

3 Method 

3.1 Seed 

A track candidate found in the pattern tracker provides the starting seed, which 
gives already rather precise estimates of the parameters x, y, t x = p x /p z and 
t y = Py/pz at the entrance of the pattern tracker, which is at z ~ 700 cm. 
We use the inverse momentum signed according to charge, Q/p, as momentum 
parameter because it is most convenient to use in the numerical treatment of 
transport within the magnetic field. Since iterations are to be avoided, it is 
essential to have already an estimate of the momentum parameter which is good 
enough to give reliable parameter derivatives. Such an estimate can be obtained 
from the visible deflection by the magnetic field assuming that the track comes 
from the target region: 

Q [x — (z — Z t arget)t x } \J 1 + t x 

P I Byd£ • (^Zffiagnet ^target) \J ^- t x 

where / B y di is the average field integral of the main bending component, z magnet 
and Ztarget are the z coordinates of magnet center and target and x and t x are 
the impact parameter and track angle in the bending plane as given at z by the 
seed. The quality of this estimate for muons and pions from the golden B decay 
is shown in fig. [3[ The estimate gives obviously reliable results in most cases. 
The relative precision is ~ 6% in the muon and m 7% in the pion case, in the 
latter the estimate is slightly diluted by the displacement of the Kg decay vertex 
which effects the tail in fig. [3]d. The initial covariance matrix diagonal element 
of the momentum parameter is set accordingly. 

3.2 Modifications to the Concurrent Track Evolution method 

After the initial seed definition from the pattern tracker, the Concurrent Track 
Evolution mechanism JU propagates the candidate as far as possible upstream 
within the magnet tracker, employing the Kalman filter technique. This is 
achieved by a succession of growth cycles, which explore the possible continu- 
ations from layer to layer. A certain number of branches is propagated concur- 
rently to explore the available paths. We summarize here briefly the steps which 
are explained in detail in ||: 
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1. Based on the position of the last hit of a candidate, generate a list of 
detector plane parts in which the next logical hit is expected, using the 
domain navigation method described in |1J. Define a search window in each 
plane, using the Kalman filter prediction from the last hit. 

2. Create a new branch for each hit found in each search window. 

3. Filter each hit and accept the new branch if the \ 2 contribution is not too 
large. 

4. At the end of each growth cycle, locate the best candidate according to a 
quality estimator function. Discard all branches whose quality estimates 
differ by more than a given maximum value, keep only the lZ m ax best can- 
didates. 

5. Repeat the steps above, until no further growth is possible. This may be 
either because the end of the tracking system is reached, or because no 
further suitable hits are found. 

6. Select the best remaining candidate, and store it if its hit count exceeds a 
certain minimum. 

Compared to the algorithm decribed in |L], the following changes and exten- 
sions have been made: 

1. The algorithm is applied in three dimensions. The same navigation strategy 
as in [|lj is used, but a domain window is immediately discarded if it is not 
intersected by the extrapolated track candidate. 

2. Hits are considered for filtering if they are within either three standard 
deviations or 2 cm of the predicted projected coordinate u = x cosa st — 
y sina si (where a st is the stereo angle). 

3. After filtering, the new candidate is accepted if the \ 2 contribution is less 
than 16. 

4. In the quality estimator Q (see |J), which is calculated from the number of 
hits and their x 2 contributions, the weight of the latter (w^) is set to 0.05. 

5. If at least one new track candidate is obtained through successful contin- 
uation, its parent track candidate is excluded from further propagation, 
This is unlike the strategy in the projection-based tracking in without 
knowledge of the orthogonal coordinate, the probability to pick up a wrong 
hit within the prediction window was higher, so that the possibility of a 
missing hit (fault) had to be considered even if a continuation hit had been 
found. 
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6. Since the magnet part of the inner tracker has only about half as many 
layers as the outer tracker, the inner tracker hits are weighted by a factor of 
two when the number of hits is calculated. A track is accepted for storage 
if its weighted number of hits, acquired in the magnet tracker is at least 5. 

The maximum number of missing hits in sequence {N™°* lts ) is kept at 2, and the 
maximum number of track candidates treated concurrently (TZ max ) remains at 5. 



3.3 Transport of track parameters within the magnetic 
field 

Compared to the field- free case discussed in where two-parameter states are 
transported linearly, the transport within the inhomogeneous field is much more 
critical because there are five track parameters involved, and both the param- 
eters and their covariance matrix must be transported. Regarding the number 
of track candidates which appear in the course of Concurrent Track Evolution, 
very efficent numerical procedures are needed in order not to exceed practical 
cpu time limitations. 

It has already been shown that the precision requirements of HERA-B track- 
ing can be met using a fifth-order Runge-Kutta method with adaptive step size 



control ||, a procedure which has also been successfully used for track fitting [[LG 
For the track propagation problem discussed in this paper, we have applied the 
following flexible strategy (described in detail in |TIj]) to keep the computational 
effort at a minimum: 

short distances: for 5z <20 cm, a parabolic expansion of the trajectory, based 
on the field vector at the starting point, has been used. 

medium distances: for 20 cm< 5z <60 cm, a classical fourth-order Runge- 
Kutta method |[2| has been applied. 



long distances: for Sz >60 cm, a fifth-order Runge-Kutta method with adap- 



tive step size control has been employed. 



Since the precision for the derivatives of the transported parameters (x, y, 
t x , t y , Q/p) with respect to the untransported ones (x, y, t x , t y , Q/p) is less 
critical than the precision of the parameters themselves, the sizes of the deriva- 
tives were studied and their influence on the track finding result was carefully 



investigated [11 . It was found that the derivatives with respect to the first four 



parameters can be approximated by 

d(x, y, t x , t y , Q/p)/dx = (1, 0, 0, 0, 0) 

d{x, y, t x , t y , Q/p)/dy = (0, 1, 0, 0, 0) 
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d{x, y, t x , t y , Q/p)/dt x = {z-z, 0, 1, 0, 0) 

d(x,y,t x ,ty,Q/p)/dt y = (0,z - z,0, 1,0) 

without loss in efficiency but with large gain in speed. The derivatives with 
respect to the momentum parameter are approximated as 

d(x,y,t x ,t y ,Q/p)/d(Q/p) = (dx/d(Q/p),0,dt x /d(Q/p),0,l) 

where dx/d(Q/p) and dt x /d(Q/p) are defined by differential equations which are 
solved together with the integration of the equations of motion by the Runge- 
Kutta method. 

These approximations lead to a derivative matrix whose elements contain 
mainly ones and zeroes. This property is used to reduce strongly the number of 
computations for the transport of the covariance matrix. 

The propagation of electrons and positrons within the magnetic field is com- 
plicated by the energy loss through bremsstrahlung. An additional radiative 
energy loss correction according to the method of Stampfer et al. [13, [?J has been 



applied, which was found to improve the magnet tracking efficiency for electrons 
from the golden B decay by 0.5% with respect to the number of reference tracks 
(sec. 



4 Performance 



4.1 Event sample 

Samples of 1219 interactions of the type pAl —> B° + X, B° —>■ J/ipKg — > 
(£ + £~)( y 7T + 'n-~), with 915 decays into the muon and 304 into the electron chan- 
nel, and 2000 unbiased inelastic pAl interactions were generated |TJJ and passed 
through the full HERA-B detector simulation |15| to investigate the performance. 
Both the event generation and the detector digitization were done in the same 
manner as described in [HI. 



4.2 Performance estimators 

As in |1[], efficiencies are calculated for reference tracks which are a track class 
defined by both the geometrical acceptance and the physics interest of the ex- 
periment. Since a track segment from the pattern tracker is needed as a starting 
point, a magnet reference track has to satisfy all criteria for a reference track in 
the pattern tracker including the momentum requirement p > 1 GeV/c. In 
addition, the track is required to pass at least four layers within the area of high 
magnetic field, defined by a circle of 140 cm radius around the magnet center in 
the bending plane. It is also demanded that the track pass at least one layer in 
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the first half of the magnet (z < 450 cm), to reject late decays where the particles 
do not traverse a sufficient amount of integrated field. Since it can happen that 
particles at the outer and inner borders of the acceptance "circumvent" some 
superlayers, tracks with such gaps are also regarded as not being covered by the 
geometrical acceptance. Finally, tracks are not considered as belonging to the 
reference set if they have lost more than 30% of their energy within the magnet. 
This requirement is directed against radiating electrons which have a reduced 
probability to be selected by the first level trigger. 

A track candidate is regarded to properly reconstruct a particle if at least 70% 
of its hits were caused by this particle. In addition, it is checked that the first 
reconstructed hit of the candidate is not further than 100 cm downstream of the 
first hit of the real particle. The efficiency is then simply the fraction of reference 
tracks which are successfully reconstructed. In order to isolate the effect of the 
magnet propagation, we define in addition a magnet pattern recognition efficiency 
which is defined relative to the track finding efficiency in the pattern tracker: 

_ ^magnet+pattern tracker 

^-magnet 

^-pattern tracker 

The efficiency e pa ttem tracker has been re-evaluated for the geometry and simulated 
data sample used here and was found to agree, within statistical errors, with the 
numbers presented in [[[]]. 

4.3 Geometrical acceptance 

In order to verify the adequateness of the reference track criteria, fig. |] shows 
the number of layers passed in the magnet, as defined in the last section, as 
a function of pseudorapidity n for muons from the golden B decay which are 
reference tracks in the pattern tracker. In the outer tracker dominated regime 
around 77 = 3, the layer count is higher than in the region governed by the 
inner tracker around n = 4.5. Obviously the cut NL ayers > 4 gives a generous 
description of the geometrical acceptance. The nominal rapidity coverage 0, 
corresponding to 10 mrad at the inner and 250 mrad (horizontal) resp. 160 mrad 
(vertical) at the outer border, is displayed as a pair of vertical lines in fig. [|. 

The fraction of particles which are selected as reference tracks can be in- 
terpreted as the geometrical acceptance of the combined system of pattern and 
magnet tracker. Table [I] summarizes these values for particles related to the 
golden B decay. Comparison of the muon acceptance of 80.2% with the one of 
the pattern tracker alone |1| shows that 92% of the muons within the pattern 
tracker acceptance have also enough hits in the magnet tracker. The electron ac- 
ceptance is smaller (66.6%) because of the cut on the radiated energy. The pion 
and K s acceptances are mainly limited by the Kg decay length which causes a 
considerable fraction of kaons to decay too late for a good momentum measure- 
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ment of the pions. The total fraction of golden B° decays with reconstructable 
condition is 34% in the muon and 29% in the electron channel. 

4.4 Efficiency 

Figure [5] shows displays of an event with seven superimposed interactions with 
focus on the magnet tracker region. The Monte Carlo tracks have been drawn by 
connecting the impact points in the sensitive detector parts, note that because of 
this technique, low-momentum particles have a poor representation. The region of 
highest magnetic field which is used for the definiton of reference tracks is shown 
as a filled circle. While several softer particles are not reconstructed because 
they do not reach the pattern tracker, the vast majority of tracks passing the full 
spectrometer is found in full length up to the entrance of the magnet tracker. 
The second display (fig. Qd) is restricted to the particles from the golden B decay 
and shows the proper recognition of these tracks. 

The pattern recognition efficiencies obtained in B interactions with on average 
four superimposed inelastic interactions are summarized in tab. |2|. The right 
column gives the relative efficiency of the magnet propagation step. For charged 
particles from the golden B decay, these efficiencies are larger than 97%. The 
tt^ appears to have a slightly lower efficiency than the leptons which may be 
attributed to the late decays of some K^s and the resulting smaller number of 
hits, and the on average smaller momentum. Charged particles above 1 GeV/c 
have a magnet propagation efficiency of ~ 94%. 

The left column shows the efficiency of the combined pattern and magnet 
tracker reconstruction. The ratio to the numbers in the right column corresponds 
within statistical errors to the numbers given in [|lj. It is interesting to note that, 
after the requirement that reference tracks should not radiate more than 30% 
of their energy within the magnet, electrons and muons have similar efficiencies. 
The relative rate of ghostsQis 2%, corresponding to an average of 0.85 ghost tracks 
per event. This is much smaller than the 4.4 ghosts per event which remained 
after pattern tracker analysis 0, which shows the ghost rejection power of the 
magnet tracker. Further suppression of ghosts is expected from the matching 
with the vertex detector system. 

If one includes the geometrical acceptance, the total fraction of reconstructed 
golden B decays is 28% in the muon and 24% in the electron channel. 

In order to probe the robustness of the algorithm, the efficiency of the magnet 
tracker propagation step has been studied as a function of the number of super- 
imposed inelastic interactions. Figure |] shows that the magnet tracking efficiency 
remains at the same high level up to ten superimposed interactions, so that the 
total efficiency is governed by the behaviour from the pattern tracker. The ghost 
rate increases in a well-controlled way. 

5 relative to the number of reference tracks 
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In order to investigate further the causes for particle loss, we have broken down 
the magnet efficiency as a function of momentum for muons and pions from the 
golden B decay (fig. [7]). For both particle species, the efficiency saturates on a 
high level above 5 GeV/c, below this value the efficiency drops, which is likely 
to be due to multiple scattering effects, or to strongly curved trajectories which 
are not entirely covered by the navigation tables. At the bottom of fig. |7|, it is 
indicated how the different parts of the momentum range contribute to the mean 
inefficiency, 1 — e Magnet- While muons and pions show the same efficiency at a 
given momentum, the pions suffer more strongly from the reduced efficiency at 
lower momentum because of their softer spectrum. 

Figure [8] shows a similar distribution as a function of the polar angle, measured 
at the track point of lowest z within the magnet tracker. The magnet track finding 
efficiency decreases visibly with increasing angle. Since, however, polar angle 
and total momentum are strongly correlated within the HERA-B kinematics, the 
effect is not independent from the one in fig. [7]. 



4.5 Parameter estimates 

The underlying reconstruction concept |IJ foresees that the track candidates 
emerging from the magnet track reconstruction should be matched with track 
candidates found in the vertex detector, or continued into the latter by direct 
propagation, which is only possible if the track parameters are well estimated. 
In principle, a fully iterative refit of the trajectory to the hits found by pattern 
recognition, which takes also additional traversed material into account, can be 
used to optimize the parameter estimate. Such a fit, the method of which is 



described in |10| , has been used in the following for comparison. From the aspect 
of cpu time consumption, it would however be preferred if a global refit can be 
postponed until the hits from all components have been collected. 
Figure |9] shows the normalized residuals 

p^REC _ %mc 



r l/2 
u ii 

where Xf- EC is the reconstructed parameter at the point of lowest z of a track 
candidate, X^ c the corresponding Monte Carlo-truth and Cu the estimate for the 
corresponding covariance matrix diagonal element, where the parameters are X^ — 
x, y, t x and t y . (If several track candidates exist for the same particle, the better 
one is chosen because it is more likely to be selected in the subsequent matching 
with the vertex detector.) The error estimate Cj/ 2 which is used to normalize the 
residual is determined by the coordinate resolution and the multiple scattering 
effects as they are estimated during magnet propagation, or through the full refit. 
Dilutions due to wrong hits or left/right assignments in the pattern recognition 
process are expected to lead to widths larger than one. All observed residual 
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distributions have shapes close to Gaussians. There are hardly any tails, as is 
also reflected in the underflow and overflow percentages which are given in the 
diagrams, which shows that the criterion used to label a track as reconstructed 



(see section |4.2|) is adequate. The fitted widths exceed the ideal value by 17- 
27%, which can be attributed to inevitable pattern recognition effects. Only in 
the case of the t x parameter the excess is smaller (4%), which is not unexpected 
because the resolution in this parameter is dominated by multiple scattering. 
One concludes that the track parameters delivered by the magnet propagation 
are good enough to be directly used by the matching step without a refit at this 
stage. 

The quality of the momentum parameter estimate influences the matching to 
the vertex detector segment only if the first point (at lowest z) is well within 
the magnetic field, because otherwise a linear extrapolation should be sufficient. 
Figure [It] shows in its upper part the normalized residual of the momentum pa- 
rameter (Q/p), which does also show only moderate tails and a width only slightly 
larger than unity. The bottom part of the same figure displays the distribution of 
the relative momentum error, Ap/p. Note that the distribution is not expected 
to be Gaussian because of the momentum spread and multiple scattering effects. 
The mean relative resolution obtained by fitting the central part of this distri- 
bution is 8.7 • 10~ 3 , which the refit improves slightly to 8.1 • 10~ 3 . To compare 
these resolutions with the technical limit, the same events were passed through 
a full iterative refit of the hits selected using Monte Carlo informationF], where a 
mean resolution of 7.1 ■ 10~ 3 was obtained. The resolutions are summarized in 
table ^. The comparison thus gives an illustration of the size of dilutions caused 
by pattern recognition effects, which turn out to be relatively moderate. Further 
improvement can be expected from 

1. a rejection of outliers, based on the y 2 contribution of each hit to the final fit, 
which can diminish the dilution from wrong hits or left/right assignments. 

2. successful matching with vertex detector segments, which should reduce the 
relative influence of wrong hit information within the magnet. 

We conclude that the method studied here gives reasonable parameter estimates 
in spite of the high track densities. 



4.6 Speed 

In the magnet case, the issue of computing time is potentially even more criti- 
cal than in the field-free case because of the costly numerical parameter trans- 
port within the magnetic field. The time consumption of the magnet part has 
been tested by running the reconstruction algorithm with and without magnet 

6 This procedure is sometimes called ideal pattern recognition 
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propagation step and computing the difference. The typical precision of such 
measurements is ~ 20%. The resulting behaviour as a function of the number 



of superimposed interactions in displayed in fig. 11. On an SGI challenge work- 



station, the computing time in case of on average four superimposed inelastic 
interactions is about 3.6 sQ The increase with the number of interactions is not 
much steeper than linear, corresponding to an almost constant time requirement 
per track. Again, the speed behaviour of the Concurrent Track Evolution algo- 
rithm turns out to be favourable for online processor farm applications, which 
was already observed in the field- free system Jl|. 



5 Summary 

The Concurrent Track Evolution algorithm has been extended and applied to the 
problem of associating detector hits with tracks in the inhomogeneous magnetic 
field of a forward B spectrometer. The method has been tested on simulated 
events with the currently implemented HERA-B geometry. The efficiency for 
propagating tracks from the golden B decay was found to be 97% and better. 
The ghost rate from the pattern tracker segments is considerably reduced by 
the magnet propagation step. The resulting track parameters are good enough 
for the vertex detector match without additional iteration or refit. The cpu 
time consumption is smaller than for the initial pattern tracker reconstruction 
step, and shows an uncritical behaviour with increasing number of superimposed 
interactions. 



Acknowledgments 

We would like to thank Thomas Lohse for contributing his tracking expertise 
in many inspiring discussions, and Siegmund Nowak for his continuous support 
regarding the HERA-B detector simulation software. One of us (A.S.) would like 
to thank DESY Zeuthen for the kind hospitality extended to him during his visit. 
This work was supported by the Bundesministerium fur Bildung, Wissenschaft, 
Forschung und Technologie under contract number 05 7 BU 351 (5). 



References 

[1] R. Mankel, A Concurrent Track Evolution Algorithm for Pattern Recognition 
in the HERA-B Main Tracking System, Nucl. Instr. and Meth. A 395 (1997) 
169. 

7 It is expected that the processors of the HERA-B reconstruction farm will be faster by at 
least a factor of 2. 



11 



[2] T. Lohse et al., An Experiment to Study CP Violation in the B System Using 
an Internal Target at the HERA Proton Ring (Proposal), DESY-PRC 94/02 
(1994). 

[3] S. Amato et al., LHCb - a Large Hadron Collider Beauty Experiment for 
Precision Measurements of CP-Violation and Rare Decays (Technical Pro- 
posal), CERN/LHCC 98-4 (1998). 

[4] R.E. Kalman, Trans. ASME, J. Basic Engineering (1960); 
R. Battin, Am. Rocket Soc. 32 (1962) 1681; 

R.E. Kalman and R.S. Bucy, Trans. ASME, J. Basic Engineering (1962); 

[5] R. Fruhwirth, Nuclear Instr. and Meth. A262 (1987) 444. 

[6] P. Billoir and S. Qian, Nuclear Instr. and Meth. A294 (1990) 219. 

[7] R. Mankel, Application of the Kalman Filter Technique in the HERA-B 
Track Reconstruction, HERA-B Note 95-239 (1995). 

[8] E. Hartouni et al., An Experiment to Study CP Violation in the B System 
Using an Internal Target at the HERA Proton Ring (Design Report), DESY- 
PRC 95/01. 

[9] T. Oest, Particle Tracing through the HERA-B Magnetic Field, HERA-B 
Note 97-165 (1997). 

T. Oest, Fast Particle Tracing through the HERA-B Magnetic Field, HERA- 
B Note 98-001 (1998). 

[10] R. Mankel, ranger - a Pattern Recognition Algorithm for the HERA-B Main 
Tracking System, Part IV: The Object-Oriented Track Fit, HERA-B Note 
98-079 (1998). 

[11] A. Spiridonov, Optimized Integration of the Equations of Motion of a Particle 
in the HERA-B Magnet, HERA-B Note 98-133 (1998). 

[12] W. Press, S. Teukolsky, W. Vetterling and B. Flannery, Numerical Recipes 
in C; Cambridge University Press (1992) 

[13] D. Stampfer, M. Regler and R. Fruhwirth, Comput. Phys. Commun. 79 
(1994) 157-164. 

[14] T.S. Sj6strand, CERN-TH-6488-92 (1992); 

B. Anderson, G. Gustafson and Hong Pi, Z. Phys. C57 (1993) 485. 

[15] S. Nowak, HBGEAN and HBRCAN, HERA-B Internal Note (1995). 



12 



Tables 



Particle 


Geometrical Acceptance 




(80.2 ± 1.3)% 
(66.6 ±1.9)% 
(66.8 ±1.1)% 


— > e + e _ 

-»• TT + TT- 


(63.4 ±1.6)% 
(45.4 ±2.9)% 
(54.8 ± 1.6)% 


5 U -»■ Jji)K% -> /i + /i-7r+7r- 
5° -»• J/^Jfg -> e+e-TT+Tr- 


(34.0 ±1.6)% 
(28.6 ±2.6)% 



Table 1: Geometrical acceptance for particles related to the golden B decay. 



Particle 


e 

Pattern ± Magnet 


e 

Magnet 




(96.8 ± 0.9)% 
(97.3 ± 0.4)% 
(93.3 ± 0.7)% 


(98.7 ± 1.1)% 
(99.6 ± 0.6)% 
(97.4 ± 0.9)% 


J/ip — > e + e~ 
X£ -> 7T+7r- 


(93.5 ±2.1)% 
(94.5 ±0.9)% 
(88.0 ±1.5)% 


(97.3 ±2.6)% 
(99.3 ± 1.3)% 
(94.8 ±1.9)% 


5° -»■ J/^s -> e+e-TT+Tr- 
5° -»• J/^Kg -> /i+/i-7r+7r- 


(83.9 ±3.9)% 
(83.0 ±2.1)% 


(92.3 ±4.9)% 
(93.6 ±2.9)% 


X ± (p > 1 GeV/c) 
Ghosts 


(88.1 ± 0.2)% 
(2.0± 0.1)% 


(94.1 ± 0.2)% 



Table 2: Pattern recognition efficiency for ordinary particles and those related to 
the golden B decay. The relative ghost rate is also shown. 
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Method 


Ap 

Gaussian fit er 


lv 

FWHM 


magnet propagation 


(8.66 ± 0.27) -1(T 3 


(20.4 ± 0.6) -10- 3 


magnet propagation 
+ refit 


(8.08 ± 0.29) -10- 3 


(19.0 ± 0.7) -lO- 3 


ideal pattern recognition 
+ fit 


(7.12 ± 0.22) -lO" 3 


(16.8 ± 0.5) -lO- 3 



Table 3: Relative momentum resolution Ap/p determined both from fitting a 
Gaussian to the central peak of the residual distribution (middle column), and 
as the full width at half maximum (right column). The resolution is given as (a) 
directly obtained from the magnet propagation procedure, (b) from a refit to the 
hits returned by this procedure and (c) from a fit to the hits selected using Monte 
Carlo information (ideal pattern recognition + fit). 
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Figure 1: Layout of the HERA-B main tracking system (top view) in front of 
the RICH. The proton beam enters from the left side. The tracking areas used 
for track finding and track propagation (see text) are indicated at the bottom. 
The shaded areas are the pole shoes of the magnet (left) and the first part of 
the RICH (right). This paper concentrates on the track propagation from the 
pattern tracker into the magnet tracker. 
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Figure 2: Magnetic field components as function of the z coordinate in the magnet 
center (a) and at a vertical displacement of 40 cm (b). The z coordinate is given 
relative to the magnet center at z magnet = 450 cm. 
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Figure 3: Correlation (a-b) between the momentum estimate (signed according 
to charge) from the apparent track deflection (see text), and the corresponding 
Monte Carlo value for muons (a) and pions (b) from the golden B decay, (c-d) 
show the distribution of the relative momentum error for the muons (c) and pions 
(d). 
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Figure 4: Number of layers passed in the magnet as a function of pseudorapidity 
r). The different levels of layer count around rj = 3 and rj = 4.5 reflect the different 
number of superlayers of the outer and inner tracking systems. The vertical lines 
correspond to the nominal acceptance region as described in the text. 
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Figure 5: (a) Display of an event with one interaction containing the golden B 
decay and six superimposed inelastic interactions, focussed on the magnet area. 
Both the Monte Carlo tracks (light grey) and the reconstructed tracks (thick dark 
lines) are show (reconstructed hit points denoted by crosses), (b) Same event, 
with the display restricted to particles from the golden B decay. 
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Figure 6: Pattern recognition efficiency in the magnet versus number of superim- 
posed inelastic interactions, for particles which traverse at least five layers in the 
central part of the magnet, and have been detected as segments in the pattern 
tracker. The ghost rate is also shown. 
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Figure 7: Magnet pattern recognition efficiency as a function of momentum for 
muons (filled circles) and pions (open circles) from the golden B decay. At the 
bottom, the contribution of each momentum bin to the mean inefficiency 1 — 
^Magnet (multiplied by 10) is shown, indicating that the smaller pion efficiency is 
a result of the softer momentum spectrum. 
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Figure 8: Magnet pattern recognition efficiency as a function of polar angle for 
muons (filled circles) and pions (open circles) from the golden B decay. At the 
bottom, the contribution of each angle bin to the mean inefficiency 1 — ^Magnet 
(multiplied by 10) is shown. 
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Figure 9: Normalized residuals (see text for explanation) for the parameters x, y, 
t x and t y for muons from the golden B decay. The solid points with error bars are 
directly obtained from the parameters delivered by the magnet propagation, the 
shaded histogram is the result of an iterative refit which interpolates traversed 
material between the hits. The result of Gaussian fits to both distributions is 
also shown, and the standard deviation a is quoted for the refitted case. The 
percentages at the arrows in the lower parts of the figures give the number of 
entries below and above five estimated standard deviations for the parameters 
delivered by the refit. 
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Figure 10: Distributions of (a) the normalized residual for the parameter Q/p, 
and (b) the relative momentum residual Ap/p for muons from the golden B decay. 
The solid points with error bars correspond to the parameter as it is delivered by 
the magnet propagation, the shaded histogram is the result of the iterative refit. 
The results of Gaussian fits to both distributions (for (b) to the central part) is 
superimposed, where the standard deviation a is quoted for the refitted case. 
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Figure 11: Computing time of the magnet propagation step per event as a func- 
tion of the number of inelastic interactions which are superimposed on the inter- 
action generating the golden B decay. 
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